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1. Black-Scholes implied volatility 
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2. Risk-adjusted call price minus Black-Scholes 
call price 
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4. Risk-adjusted gamma minus Black-Schoies 
gamma 




#include <math.h> 
#include <stdlib.h> 
#include <iostream. h> 



#define MAX_TIME_STEPS 200 

#define MAX_PRICE_STEPS 200 

#define pi 3.1415962 

#define PUT 'P' 

#define CALL 'C 

#define mxjSMMA 128.0 



#define max(a,b) (((a) > (b) 
#define min{a,b) (((a) < (b) 



(a) 
(a) 



(b) ) 
(b) ) 



double 0[]y[AX_TIME_STEPS] [MAX_PRICE_STEPS] ; // 



double 

double 

double 

double 

double 

merits 

.int SO 



Delta [MAX_TIME_STEPS] [MAX_PRICE_STEPS ] 
Gamma [MAX_TIME_STEPS] [MAX_PRICE_STEPS] 
Theta [MAX_TIME_STEPS] [MAX_PRICE_STEPS] 
Vega [MAX_TIME_STEPS] [MAX_PRICE_STEPS] ; 
M[MAX_TIME_STEPS] [MAX_PRICE_STEPS ] ; // 



array of the option price 

// array of deltas 

// array of gammas 

// array of thetas 
// array of vegas 
array of the volatility adjust 



^fbid RAPoptionEuropean ( // option's characteristics: 

T is time to maturity in years 

K is the strike price in dollars 

S is the current underlying price in dollars 

Y is the option type: Put or Call 





double T, 


// 




double 


// 




double S, 


// 




char Y, 


// 




// model param 




double sig. 


// 




double mu, 


// 




double R, 


// 




double C, 


// 




double V, 


// 




double g. 


// 




double r. 


// 




double rb, 


// 




double rd, 


// 



// computational parameters: 
int nt, // number of time steps 
int ns, // number of price steps 



double Is, 
double us 
) { 



// lower bound on underlying price S 
// upper bound on underlying price S 



// DO 
double 
eter 
double 
double 
double 
double 
prices 
double 
double 
double 
double 
double 



: initialize the computational variables 

ml = pow( 324.0 * R * C * C / pi, 1.0/3.0 ); // m is the risk param 



dt = T / nt; 
ds = ( us - Is ) / ns; 
s [MAX_PRICE_STEPS] ; 
s2 [MAX_PRICE_STEPS] ; 

t [MAX_TIME_STEPS] ; 
o[MAX_PRICE_STEPS] ; 
delta [MAX_PRICE_STEPS] 
gamma [MAX_PRICE_STEPS] 
theta [MAX_PRICE_STEPS] 



// time step 

// the price step size 

// array of the underlying stock prices 

// array of squares of the underlying stock 

// array of the time steps 

// array of the option price 

// array of deltas 

// array of deltas 

// array of deltas 



1 



^j^EPS] ; // array of vegas 



double vega [MAX_PRICEJ| 

double m[MAX_PRICE_STEPS] ; // array of the volatility adjustments 

double o_old[MAX_PRICE_STEPS] ; // array of old option prices 

int it; // an auxiliary index variable 
int is; //an auxiliary index variable 
double ds2 = ds * ds; 

double a [MAX_PRICE_STEPS] ; // array of matrix coefficients 
double b [MAX_PRICE_STEPS] ; // array of matrix coefficients 
double c [MAX_PRICE_STEPS] ; // array of matrix coefficients 
double d[MAX_PRICE_STEPS] ; // array of matrix coefficients 
double auxl, aux2 , aux3 , aux4, aux5, aux6; 

// DO: Make sure that the coefficients are within their meaningfull range 

if ( nt > MAX_TIME_STEPS | | nt < 1 ) nt = MAX_TIME_STEPS ; 

if ( ns > MAX_PRICE_STEPS j | ns < 1 ) ns = MAX_PRICE_STEPS ; 

if( sig < 0.0 ) sig = - sig; // volatility cannot be negative 

if( R < 0.0 ) R = -R; // variance risk premium cannot be negative 

if ( g < 0 . 0 ) g = -g; // gap risk premium cannot be negative 

if( rb < 0.0 ) rb = -rb; // borrowing fees cannot be nagative 

if ( Is < 0 . 0 ) Is = 0 . 0 ; // underlying price range cannot be negati 

.^e 

y if ( us < 0 . 0 ) us = -us; // underlying price range cannot be negati 

V = 0.0; // set the volatility risk premium to zero in this demonstration 
{procedure 

'"J // DO: Set the borrowing rate and dividend rate for puts to zero 
ffl if ( Y == PUT ) { 
M rb = 0.0; 

L ^ 

£ // DO: Initialize the arrays of the time 

□ for( it = 0; it <= nt; it++ ){ 

U t[it] = it * dt; // the maturity time T is t[nt], the current time 0 is 

// DO: Initialize the arrays of the price 
is = 0; 

s [0] = S - is * ds; 
while( s[0] >= Is ) { 

s[0] = s[0] - ds; 

is++; 

} 

SO = is - 1; 
s [ 0 ] = s [ 0 ] + ds ; 
for( is = 1; is <= ns; is++ ){ 
s[is] = s[is-l] + ds; 

} 

for( is = 0; is <= ns; is++ ){ 

s2[is] = s[is] * s[is]; // array of the price squares 

} 

// DO; initialize the finite difference scheme by the derivative's values 
// at maturity time by the payoff function 

if ( Y == CALL ) { 

for( is = 0; is <= ns; is++ ){ 

o[is] = max( 0.0, s[is] - K ); // call option payodeltasff 
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if( s[ 
delta [is 
gamma [ i s 

} else 
delta [is 
gamma [ i s 

} else 
delta [is 
gamma [is 

} 

theta [ 
vega [i 

} 

else if( 
f or ( is 
o [is] 
if( s[ 
delta [is 
gamma [ i s 
} else 
delta [is 
gamma [is 

} else 
delta [is 
gamma [ i s 

} 

theta [ 
vega [i 



is] < K ) { 

] =0.0; // delta zero below the strilce 
] =0.0; // gamma zero 

if ( s[is] == K ) { 
] =0.5; //at the strike delta 0.5 for numerical stability 
] = MAX_GAMMA; // gamma is the maximum allowed by the marlcet 

// which is inverse of the minimum ticlc size 

if ( s[is] > K ) { 
] =1.0; // delta one above the strilce 
] =0.0; // gamma zero 

is] =0.0; //no time decay left 

s] = 0.0; // no volatility sensitivity 

Y == PUT ) { 
= 0; is <= ns; is++ ) { 

= max( 0.0, K - s[is] ); // put option payoff 
is] < K ){ 

] = -1.0; // delta zero below the strilce 
] =0.0; // gamma zero 

if ( s[is] == K ) { 
] =0.5; //at the strilce delta 0.5 for numerical stability 
] = MAX_GAMMA; // gamma is the maximum allowed by the market 

// which is inverse of the minimum tick size 

if { s[is] > K ) { 
] =0.0; // delta one above the strike 
] =0.0; // gamma zero 

is] =0.0; //no time decay left 

s] = 0.0; //no volatility sensitivity 



} 

for( is = 0; is <= ns; is++ ){ 
if ( gamma [is] > 0.0 ) { 

m[is] = max( 0.0, 1.0 + g - min ( ml 

2 .0 * R * delta[is] 

} else { 

m[is] = 1.0 + g; 



* pow{ gamma [is], 1.0/3.0 
* delta[is] / gamma[is] ) 



} 



} 



for( is - 0; is <= ns; is++ ){ 



} 



0[nt] [is] 
Delta [nt] [is] 
Gamma [nt ] [is] 
Theta[nt] [is] 
Vega[nt] [is] 



o[is] ; 
delta [is] 
gamma [is] 
theta [is] 
vega [is] ; 



// DO: remember the last computed values 
for( is = 0; is <= ns; is++ ){ 
o_old [is] = o [is] ; 

} 

// DO: Setting: starting time is maturity time 
it = nt; 



// DO: set the auxiliary variables for faster computational speed 
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dt / ( ds * ds ) ; 



auxl = 0.25 * sig 
aux2 = 0.25 * ( r - rb - rd ) * dt / ds; 
aux3 = 0.5 * r * dt; 

// DO: in the loop compute all the optionis prices 
while ( it >= 0 ){ // INV: new time >= current time 

it--; // Meaning: Time := Time - Time Step 

// DO: compute the solution on the new time level by solving 

// the RAP equation with the semi-implicit finite difference 

// Crank-Nicolson scheme method 

// DO: Initialize the triangular matrix coefficients and the right hand 

side 

f or ( is = 1; is < ns; is++ ){ 
// a,b,c are on the diagonal 
aux4 = auxl * s2 [is] * m[is] ; 
aux5 = aux2 * s[is]; 

a [is] = aux4 - aux5; 
b[is] = -1.0 - 4 * aux4 - aux3 ; 
c[is] = aux4 + aux5; 
// d is the right hand side 

d[is] = o_old[is-l] * ( - aux4 - aux5 ) + o_old[is] * ( -1 + 4 * aux4 
rH aux3 ) 

W + o_old[is + l] * ( - aux4 + aux5 ) + v * vega[is]; 

} 

//DO: Initialize the boundary conditions: 
if ( y == CALL ) { 
O // Call boundary conditions 
£ b[0] = 1.0. 
c[0] = 0.0 
d[0] = 0.0 

a[ns] = 0.0; 
b[ns] = 1.0; 

d[ns] = ( s[ns] - K ) * exp ( ( r - rd - rb ) * ( T - t[it] )); 
} else if ( Y PUT ) { 

// Put boundary conditions 
b[0] = 1.0; 
c[0] = 0.0; 

d[0] = ( K - s[ns] ) * exp( - r * ( T - t[it] )); 

a[ns] = 0.0 
b[ns] = 1.0 
d[ns] = 0.0 

} 

// DO: Solve the matrix for the new options values 
// DO: run the forward elimination first 
for( is = 1; is < ns; is++ ){ 
aux6 = - a[is] / b[is-l]; 
b[is] += c[is-l] * aux6; 
d[is] += d[is-l] * aux6; 

} 

// DO: run baclcwards elimination to get the options values 
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o [ns] = d[ns] ; 
for( is = ns-1; is > 0; is-- ){ 

o[is] = ( d[is] - c[is] * o[is+l] ) / b[is] ; 

} 

o[0] = d[0] ; 

// DO: compute the sensitivities 
for( is = 1; is < ns; is++ ){ 

o[is] = max( 0.0, o[is] ); // option value should not be negative 

delta[is] = ( o[is+l] - o[is-l] ) / ( ds + ds ) ; 

gamma [is] - min( MAX_GAMiy[A, 

max( 0.0, ( o[is+l] - 2 * o[is] + o[is-l] ) / ds2 ) ); 

theta[is] = ( o__old[is] - o[is] ) / dt; 

vega[is] =0.0; //in this software vega need not be computed 

// because we assume v zero 

i f ( gamma [ i s ] > 0 . 0 ) { 
m[is] = max( 0.0, 1.0 + g - min ( ml * pow( gamma [is] , 1.0/3.0 ), 

2.0 * R * delta[is] * delta[is] / gamma[is] ) ); 

} else { 
m[is] = 1.0 + g; 
} 

// DO: save the previous step for the options 
o_old [is] = o [is] ; 

} 

// DO: Save the computed values into the arrays 
for( is = 0; is <= ns; is++ ){ 

0[it] [is] = o[is] ; 

Delta[it] [is] = delta[is]; 

Gamma [it ] [is] = gammafis]; 

Theta[it] [is] = theta[is]; 

Vega [it] [is] = vega [is]; 

M[it] [is] = m[is] ; 

} 

} 

// INV: the values have been computed up to the current time 

// Return of the computed derivative's values and sensitivities 

// is done throughthe globally defined arrays 0, Delta, Gamma, Theta, Vega 
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F.c^u/< 1- ^ Options 

No mystery 

behind 
ttie smile 

Milan Kratka presents a risk-ad|usted pricing model that 
matches the market prices of options more closely 



Classical option theory derives option values from the 
so-called arbitrage replication process: investors who 
buy options and hedge them continuously and with- 
out cost are expected to generate a risk-free return 
equal to that from Treasury bonds. These idealistic 
assumptions imply that the option model is not con- 
sfeijnt with real markets. In real life, every investor rehedges at discrete 
tihae intervals, pays trading costs, is exposed to potential price gap and 
ve^tility risks, and experiences positive portfolio variance. 

fjhe pricing model we present here helps close this gap between theory 
affS" practice. Risk-adjusted pricing incorporates all the real-life factors into 
assimple but comprehensive model that is consistent with the market. It 
ef^nces the old "no-arbitrage" paradigm with the more realistic notion 
ofSalue-at-risk. We will try to explain the main ideas behind this break- 
tHgugh methodology. 

Siandard assumptions 

Sv^ose that an investor prices an option 0 on a stock S that pays zero 
dj^ends and follows the typical, constant-volatility lognormal process: 

" dS = ^iSdt + aSdl 

where \i denotes the drift (or the trend) and o denotes the volatility. Let 
us assume we buy N options and short-sell 8N stocks as a hedge. We plan 
to rehedge this portfolio by trading in stock. Let us denote the rehedging 
time interval by A^. This should reflea the trading strategy chosen by the 
investor and, in general, will vary with the time to expiry, underiying price 
and other factors. Since no portfolio can be hedged continuously or with- 
out cost, the investor needs to specify a rehedging rule. For example, we 
might simply want to hedge in the most advantageous way. So we ask 
our quants to figure out how to do this. They may come up with strate- 
gies such as: 

□ rehedge on sunny days only; 

□ rehedge every hour; 

□ rehedge with every 1% change in the stock price; 

□ rehedge with every 0.05 change in delta; or 

□ leave rehedging strategy to the best trader. 

What would be the optimal strategy? Rehedging too frequendy accu- 
mulates high trading costs. Hedging too infrequently results in high port- 
folio variance. The risk-adjusted pricing methodology offers a natural 
solution to the problem of specifying an optimal rehedging strategy. 

Risk-actjusted pricing 

We have now set up a portfolio that will be rehedged in a time interval Aj. 
Over this period, we expect it to generate a return equal to that from risk- 
free securities, but it should also pay for the trading costs and compensate 
us for having positive portfolio variance and for taking volatility exposure 
risk. Though the transactions will occur at discrete times, we can locally 



estimate the rate of trading costs by spreading them over the time between 
the rehedging points, as is the customary practice with the dividend rate. 
We will estimate the rate of the portfolio variance similarly. 

Transaction fees, trading strategies and risk aversion differ from one 
investor to another. An option may seem to be undervalued for one in- 
vestor and overpriced for another. The market price of a liquid securi- 
ty reflects the prevailing expectations among all active investors. Who 
would not like to have a pricing model which is both in line with the 
market and reflects investor-specific conditions and expectations? The 
latter is not quite possible without the former. Risk-adjusted pricing al- 
lows us to model both. 

Let us denote the risk-free rate by r, xhe slock borrowing rate by fg, the 
expected rate of trading costs by r^Q, the expected rate of premium for tak- 
ing positive variance risk by fy^^ and the expeaed rate of premium for 
volatility exposure risk by r^^. The average change E[dP] of the portfolio 
value in an infinitesimal time increment dt should then satisfy: 

E[dP] = E[dO]-E[5dS] 

= [dtO + la^S^assOjdt + E[(asO - 5)dS] 

= r(0- 6S)dt + reSSdt + (rvAR + rjc + ro)dt 

Rehedging to the market delta d^O would make the expected value 
E [(dgO - 6)dS] equal zero, and the portfolio would have zero variance over 
the time dt. On the other hand, rehedging consistently to a different delta 
5 would make the expectation term equal to OgO - 6)LiSdt, and it would 
increase the expected portfolio variance by OsO-5)Vs2dt. For investors 
with strong trend expectations, it is optimal to trade to a trend-adjusted 
delta that differs from the market delta. 

in general, it is difficult to estimate the future trend term \i from the 
market. It is a subjective issue. If the trend was known and had significant 
value, the underiying price would have already been adjusted according- 
ly. There is much uncertainty about trends. Maybe that is why they are 
often referred to by the more moderate term "drifts". Each time we buy a 
stock, somebody sells it. Similarly, some investors might underhedge, while 
others might overhedge. For the market modelling we could simply as- 
sume that most traders hedge to the market delta. For investors with a 
strong opinion about the u^end, or even for a market reflecting these opin- 
ions, rehedging to a trend-adjusted delta should be more advantageous, 
since it should increase expected return. 

Poitfolio variance risk premium 

How much, in premium, should investors be compensated for positive 
portfolio variance risk? It is safe to say that the higher the variance, the 
higher the premium. As the variance of a normal process increases linearly 
in time, it may seem natural to expect the premium to follow the portfo- 
lio variance lineariy. This is similar to portfolio performance rating using 
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rvAR - ' 



Sharpe Ratios. If we denote the variance risk coefficient by R, the rate of 
the portfolio variance risk premium is: 

R • var(AP) 

You can also consider R as a risk aversion coefficient. The variance of 
the portfolio between two rehedging points can be estimated as: 

(AP-E[AP]f' 



var(AP) = E 



(asO(t + At)-5(t)) 4 
(ASf 



r2(AS)S(asO(t)-5(t)f(Asf] 



where AS denotes the change in the underlying price between rehedging 
times t and t+A^^ and T = d^^O denotes market gamma, as is usual. The 
rate of the portfolio variance risk premium is then: 

rvAR = SRr^a'^S^At + R(asO - dfo^S^ 

^^Jotably, the longer the time between rehedging A^, the higher the rate 
oWie variance risk premium r^^p. 

\lior a risk-averse investor with a strong opinion about the trend, the op- 
tii^jl delta to rehedge to is the one that maximises the risk-adjusted, trend- 
related rate: 

(dsO - S)(^l - r + Tb) S - R(asO - 6fo^S^ 



i'^hat delta is: 



8 = asO- 



:;=^his trend-adjusted delta makes the risk-adjusted, trend-related rate 
(yn rB)^/4Ra'^. A moderate up-trend would cause option holders to 
b(|Iiedged with a lower delta than the market delta. Trends seem to have 
a'Telativeiy high impact on options that are not ai-the-money. Whatev- 
er the direction and strength of future trends, their presence and mar- 
ket anticipation will help to enlarge the effect of higher implied volatility 
at the tails. 

Tiding costs estimator 

An investor's trading costs typically depend on the number of contracts 
traded, type of order, size of bid-ask spread, bid and ask volumes, ex- 
change fees and a transaction fees structure negotiated with a broker. Let 
us denote an aggregate trading cost per contract, when the size of a trade 
is n contracts, by C(n). It is typically a non-increasing function of a trade 
size, C'(n) < 0. The larger the trading size, the higher the total commissions 
but the lower the total cost per traded share. It does not need to be posi- 
tive, as heavy buying may move the market above the average cost of the 
trade, causing an immediate paper gain. For market modelling, we may 
assume that the prevailing aggregate trading cost per contract is constant 
C(n) = C. Investors should specify their own trading costs structure. It will 
affect investors' expected option theoretical values and their choice of re- 
hedging strategy. Trading costs alone would have a relatively small price 
effect but combined with discrete hedging, positive variance risk and volatil- 
ity exposure premiums can have a significant impact. 
The expected rate of trading costs: 

c(n|a5|)|a5| 

rTC= ^ \ ^ ' 

reflects that we will trade N I A5 I shares and spread out the cost over a time 
interv^al A^. The average amount of shares traded per option in each delta 
rehedge adjustment can be estimated using Ito's lemma: 



e[ I A5 1 ] = e[ |at5At + asSAS + |ass5(Asf I 

^aS|as6|E[|d7|] 
= aS|r|^V^ 



If the time between rehedging was long enough, the expected value 
could be calculated more precisely, and would marginally depend on the 
trend \i. The rate of trading costs can finally be expressed as: 

rTC=c(N^aS|r|VK7)^aS|r|-A. 

Notably, the longer the time between rehedging A^, the lower the trad- 
ing costs rate (jq. 

The higher the volatility sensitivity of the portfolio (which equals the 
option's vega), the higher the volatility risk we are exposed to. How should 
we model the volatility risk exposure? It will depend on our anticipation 
of volatility change, which is a subject for a separate project, a volatility 
model. We can model the volatility risk premium by = vla^^OL sa\-. re- 
flecting our portfolio sensitivity to any change in underlying volar ilii\\ The 
higher the vega, the higher the premium. The specific form of tiie volatil- 
ity risk premium coefficient should come from our volatility model . It should 
reflect the probability of the volatility change as well as the size of its prob- 
able change. 

Optimal hedging strategies 

If we believe in the efficiency of financial markets, we might suppose that 
market option prices reflect the market value of risk, average prevailing trad- 
ing costs and competitive rehedging strategies. There is never a guarantee 
that the free market will value options in any particular way. But we can as- 
sume it does, build a pricing model and use it to verify our presumption. 

Given investor-specific trading costs C(n) and an investor-specific vari- 
ance risk coefficient R, the optimal rehedging strategy is identified by min- 
imising the sum of r^^^p and Tjq, the variance premium and the trading cost 
terms: 

rvAR + rjc = SRa^S^r^At + C (^^Na S | r | V^]^aS | F | 

where the rehedging time A^ varies in a reasonable range. This may span 
a fraction of an hour to a fraction of a year. The gamma, T, should be equal 
to the market gamma B^O, The simple task of finding the minimum can be 
done approximately or even precisely. There is no need to write elaborate 
papers about what an optimal hedging strategy in the presence of trading 
costs ought to be. If the trading costs per share, C, do not depend on trade 
size, the optimal rehedging time would be: 

At =- 



2s2p2/3 



where we have denoted the constant: 



k = (C / SRyfOn) 



1/3 



This translates to a very simple strategy of rehedging with every stock 
change of: 

AS = ±aS7A7 = ±kr~^^^ 

This is probably close to a rehedging rule you have already used. If not, 
feel free to upgrade now! 

In line with our intuition, we shall trade more often with higher risk 
aversion R, lower trading costs C or higher gamma f. In the limiting case 
of zero transaction costs, we would trade continuously, as assumed by the 
Black-Scholes model. 

The assumption of competitive hedging in the market will then yield 
the sum of the variance premium and the trading costs terms ,to be: 



I'VAR = ^JC ~J 



ma 



2s2ppl/3 
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1. Black-Scholes implied volatility 2. Risk-adjusted call price minus Biack-Schoies 

call price 




3. Risk-adjusted delta minus Black-Scholes delta 4. Risk-adjusted gamma minus Black-Scholes 

gamma 




^:where we have denoted the constant m = {324RC2/rt)^. This was die last 
igece needed to finalise a surprisingly simple risk-adjusted pricing equation: 
atO + |a2s2(l - jp ^ r(0 - 5S) + TbSS 

It can be made more complete by including the drift term -{^t- r -I- rg)2/4Ro2 
and the volatility risk premium term that in our simplified form was v 1 9^jO I . 
It can be solved numerically backwards in time by finite difference scheme 
methods. The three-dimensional sample plots above (figures 1-4) have been 
produced by choosing an underlying volatility of 30%, a risk parameter of 
0.9 and an interest rate of 6%. For the purposes of this paper, we have left 
out the technical details of the model implementation.' We will just note 
that the rate of the portfolio variance should not be greater than the vari- 
ance S^a^gZ of an unhedged option position, which caps the variance risk 
premium that would have otherwise spiked up at the strike price at the op- 
tion's maturity. 

Volatility smile 

There are just two market parameters in die simple risk-adjusted option 
pricing equation - the anticipated volatility of the underlying price, o, and 
the so-called risk parameter m - although we would consider the volatil- 
ity risk premium coefficient v, trend parameter and risk aversion coeffi- 
cient R for a more detailed option model. The volatility risk premium 
coefficient measures our exposure to the risk that the market could change 
its mind about the anticipated volatility. Its effect on option pricing should 
diminish with time to maturity, and thus it would not explain sharpening 
of the so-called volatility smile. The risk parameter is the one responsible 
for its sharpening, as its effect on option pricing magnifies with gamma, 

' Fee/ free to contact the author if you have any questions 



as if the options with higher gamma were priced with lower volatility. The 
trendiness of the market \i helps to sharpen the smile at the tails. 

In any case, it is quite nonsensical to assume a different underiying 
volatility for different options. There can be only one anticipated under- 
lying volatility, no matter which options we look at. Even if we could in- 
corporate fat tails, price gaps, complex supply-demand relations or even 
autocon-elation, there would still be only one underiying process widi only 
one volatility. The notion of implied volatilities points out how far away 
the market pricing is from an idealised risk-neutral Black-Scholes worid. 

So, why do we observe volatility smiles? Because there are costs and 
risks associated with options trading, as we have seen and become capa- 
ble of modelling. There is no mystery behind die smile any more. The risk- 
adjusted pricing methodology may not replicate die market pricing 
perfecdy, but it does explain quite a significant part of its stable dynam- 
ics. It provides an edge in understanding option valuation and optimal 
hedging strategies. It shows how to incorporate various market risks quan- 
titatively into an option pricing equation. It will even save resources that 
would- otherwise be spent on artificial implied volatility surfaces model- 
ling. Finally, it is simple and easy to put into production. 

We have chosen to model market risks in a simplified way to demon- 
strate the natural power of the risk-adjusted pricing mediodology. It is just 
a first Step in pointing out the new dimensions possible in die financial 
modelling of a not quite risk-neutral worid. ■ 

Milan Kratka is senior partner at the Chicago-based Dragon Financial 
Group and also lectures at the University of Chicago. He would tike 
to thanic his wife ICat'lca for her constant support and Tim Weithers 
from SBC Warburg Dillon Read for his valuable comments 

e-mail: milan @math. uchicago. edu 
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